REAL- VALUED ALGEBRO-GEOMETRIC SOLUTIONS 
OF THE CAMASSA HOLM HIERARCHY 
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Abstract. We provide a treatment of real-valued, smooth, and bounded 
algcbro-geometric solutions of the Camassa-Holm (CH) hierarchy and describe 
the associated isospcctral torus. We also discuss real- valued algebro-geometric 
solutions with a cusp behavior. 



1. Introduction 

Recently, we provided a detailed treatment of the Camassa-Holm (CH) hierarchy 
with special emphasis on its algebro-geometric solutions. The first nonlinear partial 
differential equation of this hierarchy, the Camassa-Holm equation, also known as 
the dispersive shallow water equation [5] is given by 

4u t - u xxt - 2uu xxx - 4u x u xx + 2Auu x = 0, (x, t) e M 2 (1.1) 

(choosing a convenient scaling of x,t). For various aspects of local and global 
existence, and uniqueness of solutions of (1.1), wave breaking phenomena, soliton- 
type solutions ("peakons"), complete integrability aspects such as infinitely many 
conservation laws, (bi-)Hamiltonian formalism, Backlund transformations, infinite 
dimensional symmetry groups, etc., we refer to the literature provided in [12] (see 
also [13, Ch. 5]). The case of spatially periodic solutions, the corresponding inverse 
spectral problem, isospectral classes of solutions, and quasi-periodicity of solutions 
with respect to time are discussed in [7]— [11] . Moreover, algebro-geometric solutions 
of (1.1) and their properties are studied in [l]-[3], and [12] (see also [13, Ch. 5]). 

In Section 2 we recall the basic polynomial recursion formalism that defines the 
CH hierarchy using a zero-curvature approach. Section 3 recalls the stationary 
CH hierarchy and the associated algebro-geometric formalism. Section 4 provides a 
brief summary of self-adjoint canonical systems as needed in this paper, and Section 
5 finally discusses the principal result of this paper, the class of real- valued, smooth, 
and bounded algebro-geometric solutions of the CH hierarchy and the associated 
isospectral torus. 

This paper should be viewed as a companion to our paper [12] (see also [13, Ch. 
5]) on the CH hierarchy and we refer to it for background material and pertinent 
references on the subject. 



Date: August 14, 2002. 

1991 Mathematics Subject Classification. Primary 35Q53, 58F07; Secondary 35Q51. 

Key words and phrases. Camassa-Holm hierarchy, real-valued algebro-geometric solutions. 

Research supported in part by the Research Council of Norway. 

1 



2 



F. GESZTESY AND H. HOLDEN 



2. The CH hierarchy, recursion relations, and hyperelliptic curves 

In this section we review the basic construction of the Camassa-Holm hierarchy 
using a zero-curvature approach following [12] (see also [13, Ch. 5]). 

Throughout this section we will suppose the following hypothesis (No = NU{0}). 

Hypothesis 2.1. In the stationary case we assume that 

d m u 

ueC°°(R), — e J L°°(R), togN . (2.1) 
In the time- dependent case (cf. (2.28)-(2.35)) we suppose 

u(-,t) € C°° (M), g-^(-,t) G L°°(R), to G No, i G R, 

ufi,-),^,.)?^!),^!. (2.2) 

We start by formulating the basic polynomial setup. One defines {fe}eeN re- 
cursively by 

/o = l, Ax = -2a(2(4u-u xa: )/£_i, x + (4M x -u xa:a: )/^-i), feM, (2.3) 
where Q is given by 

0: L°°(M) -► L°°(R), (0v)(a;) = \ ( dy e" 2 !^^^), i£l, »6 L°°(R). 

(2.4) 

One observes that £/ is the resolvent of minus the one-dimensional Laplacian at 
energy parameter equal to —4, that is, 

d 2 \ _1 

+ 4 . (2.5) 



dx 2 

The first coefficient reads 

h = -2u + c u (2.6) 

where c\ is an integration constant. Subsequent coefficients are non-local with 
respect to u. At each level a new integration constant, denoted by eg, is introduced. 
Moreover, we introduce coefficients {ge}eeN an d {he}eeN by 

9t = ft + ^fe,x, ht = (4u-u xx )fe- ge+i, x , £eN - (2.7) 
Explicitly, one computes 

/o = l, /i = -2u + ci, f 2 = 2u 2 + 2G(u 2 x + 8u 2 ) + Cl (-2u) + c 2 , 

.90 = 1, c/i = -2u-u x + a, (2.8) 

.92 = 2u 2 + 2uu x + 2Q(u 2 x + u x u xx + 8uu x + 8u 2 ) + c\(-2u - u x ) + c 2 , 

h = 4u + 2u x , 

hi = -2u 2 x - Auu x - 8u 2 

- 2Q(u x u xxx + u 2 xx + 2u x u xx + 8uu xx + 8u 2 x + I6uu x ) 
+ ci(4u + 2u x ), etc. 
Given Hypothesis 2.1, one introduces the 2x2 matrix U by 

^ = G-w - u xx{x)) l)' x ^ ^ 
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and for each n £ No the following 2x2 matrix V n by 

Vn{z, X ) Glfc))' -£n , Z £C\ { 0}, X £M, (2.10) 

assuming F n , G n , and H n to be polynomials of degree n with respect to z and C°° 
in x. Postulating the zero-curvature condition 

-v n ) + [U(z,x),V n (z,x)} = 0, (2.11) 

one finds 

F n , x (z,x) = 2G n (z,x) - 2F n (z,x), (2.12) 

zG n , x (z,x) = (Au(x) -u xx (x))F n (z,x) -H n (z,x), (2.13) 

H niX (z,x) = 2H n (z,x) - 2(4u(x) - u xx {x))G n (z, x). (2.14) 
From (2.12)-(2.14) one infers that 

^det(V n (z,x)) = -~(zG n (z,x) 2 + F n (z,x)H n (z,x)^ = 0, (2.15) 

and hence 

z 2 G n (z,x) 2 + zF n (z,x)H n (z,x) = R 2n+2 (z), (2.16) 
where the polynomial R 2n +2 of degree 2n + 2 is x- independent, 

2n+l 

R2n+2(z) = [] ( z £ o, ^1, ... , ^2„ G C, E 2n+1 = 0. (2.17) 

m=0 

Next one makes the ansatz that F n , H n , and G n are polynomials of degree n, 
related to the coefficients ft, he, and ge by 

n n n 

F n (z,x) =^2fn-e(x)z e , G n (z, x) = ^ 9n-i{x)z e , H n (z, x) = h n _ e (x)z e . 
1=0 e=o e=o 

(2.18) 

Insertion of (2.18) into (2. 12)- (2. 14) then yields the recursion relations (2.3), (2.4) 
and (2.7) for fa and ge for I = 0, . . . ,n. For fixed n £ N we obtain the recursion 
(2.7) for hi for I = 0, . . . , n - 1 and 

h n = (4m - u xx )f n . (2.19) 

(When n = one directly gets h = {Au — u xx ).) Moreover, taking z = in (2.16) 
yields 

2n 

f n (x)h n (x) = - 11 E m . (2.20) 

m=0 

In addition, one finds 

h n , x (x) - 2h„(x) + 2(Au(x) - u xx (x))g„(x) = 0, n £ N . (2.21) 

Using the relations (2.19) and (2.7) permits one to write (2.21) as 

s-CH„(u) = (u xxx - Au x )f n - 2{Au - u xx )f n ^ x = 0, n £ N . (2.22) 

Varying n £ No in (2.22) then defines the stationary CH hierarchy. We record the 
first few equations explicitly, 

s-CHo(w) = u xxx - 4u x = 0, 
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s-CHi(w) = -2uu xxx - Au x u xx + 24uu x + c\{u xxx - Au x ) = 0, (2.23) 
s-CH 2 (m) = 2u 2 u xxx - 8uu x u xx - 40u 2 u x + 2(u xxx - <iu x )G(u 2 x + 8u 2 ) 
- 8(4u - u xx )Q(u x u xx + 8uu x ) 

+ ci(-2uu xxx - iu x u xx + 24uu x ) + c 2 (u xxx - Au x ) = 0, etc. 

By definition, the set of solutions of (2.22), with n ranging in No, represents 
the class of algebro-geometric CH solutions. If u satisfies one of the stationary 
CH equations in (2.22) for a particular value of n, then it satisfies infinitely many 
such equations of order higher than n for certain choices of integration constants 
C£. At times it will be convenient to abbreviate (algebro-geometric) stationary CH 
solutions u simply as CH potentials. 

Using equations (2.12)-(2.14) one can also derive individual differential equations 
for F n and H n . Focusing on F n only, one obtains 

(z, x) - 4(z 1 (4u(x) - u xx {x)) + l)F ntX (z, x) 
- 2z-\4u x (x) - u xxx {x))F n {z, x) = (2.24) 

and 

- (z 2 /2)F n , xx (z, x)F n (z 7 x) + (z 2 /4)F„^(z, xf + z 2 F n (z, xf 

+ z(4u(x) - u xx (x))F n (z, x) 2 = R 2n +2{z). (2.25) 

Equation (2.25) leads to an explicit determination of the integration constants 
ci,...,c„ in the stationary CH equations (2.22) in terms of the zeros Eq = 0, 
Ei, . . . , E 2 n+i of the associated polynomial R 2n +2 in (2.17). In fact, one can prove 

ce = c t (E), l = 0,...,n, (2.26) 



where 

co{K) 



c k (E) = 



(2j 1 )!...(2j 2 „ +1 V 



, : ; ; _ :0 2 2fc (ji!) 2 • • • (j 2n+ i!) 2 (2j 1 - 1) • • • (2j 2n+1 1) 

JlH h72n+l=k 

xE?---Et?l, keN. (2.27) 

Next, we turn to the time-dependent CH hierarchy. Introducing a deformation 
parameter t n G R into u (i.e., replacing u(x) by u(x,t n j), the definitions (2.9), 
(2.10), and (2.18) of U, V n , and F n , G n , and H n , respectively, still apply. The 
corresponding zero-curvature relation reads 

U tn {z,x,t n ) - V ntX {z,x,t n ) + [U(z,x,t n ),V n (z,x,t n )] = 0, n e N , (2.28) 

which results in the following set of time-dependent equations 

4u tn {x,t n ) - u xxtn {x,t„) - H n>x (z,x,t n ) + 2H n (z,x,t n ) 

- 2(4u(x, t n ) - u xx (x, t n ))G n {z, x, t n ) = 0, (2.29) 

F n ,x{zi x, t n ^j — 2G n (z, x, i n ) 2_F n (z, x : t n ), (2.30) 

zG njX (z,x,t n ) = (4u(x,t n ) - u xx (x,t n ))F n (z,x,t n ) - H n (z,x,t n ). (2.31) 

Inserting the polynomial expressions for F n , H n , and G n into (2.30) and (2.31), 
respectively, first yields recursion relations (2.3) and (2.7) for ft and ge for I = 
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0, ...,n. For fixed n £ N we obtain from (2.29) the recursion for he for t — 
0, . . . , n — 1 and 

h n = (Au-u xx )f n . (2.32) 

(When n = one directly gets ho = (4u — u xx ).) In addition, one finds 

4u tn (x,t n ) - u xxtn (x,t n ) - h n>x (x,t n ) + 2h n {x,t n ) 

-2(4u(x,t n ) - u xx (x,t n ))g n (x,t n ) = 0, n € N . (2.33) 

Using relations (2.19) and (2.32) permits one to write (2.33) as 

CH n (u) = 4u tn - u xxtn + (u xxx - 4u x )f n - 2(4u - u xx )f n , x = 0, n £ N . 

(2.34) 

Varying n £ No in (2.34) then defines the time-dependent CH hierarchy. We record 
the first few equations explicitly, 

CH (u) = Au tQ - u xxta + u xxx - Au x = 0, 

CHi(u) = 4u tl - u xxtl - 2uu xxx - Au x u xx + 2Auu x + c x {u xxx - Au x ) = 0, 

CH 2 (u) = 4u t2 - u xxt2 + 2u 2 u xxx - 8uu x u xx - 40u 2 u x (2.35) 

+ 2(u xxx - 4u x )Q(ul + 8u 2 ) - 8(4u - u xx )Q{u x u xx + 8uu x ) 

+ c 1 (-2uu xxx - 4u x u xx + 2Auu x ) + ci{u xxx - Au x ) = 0, etc. 

Up to an inessential scaling of the (x,t{) variables, CHi(u) = with C\ = repre- 
sents the Camassa-Holm equation as discussed in [5] . 

3. The algebro-geometric CH formalism 

This section is devoted to a quick review of the stationary CH hierarchies and 
the corresponding algebro-geometric formalism as derived in [12] (cf. also [13, Ch. 
5]). 

We start with the stationary hierarchy and suppose 

d m ii 

u£C oc (R), — £L°°(R), me No, (3.1) 

and assume (2.3), (2.7), (2.9)-(2.11), (2.16)-(2.22), keeping n £ N fixed. 
Recalling (2.17), 

2n+l 

R2n+2{z) — W_{z — E m ), Eo, Ei, . . . , Em £ C, E 2n +i = 0, (3.2) 

m=0 

we introduce the (possibly singular) hyperelliptic curve K, n of arithmetic genus n 
defined by 

K n : T n {z, y) = y 2 - R2n+2(z) = 0. (3.3) 

In the following we will occasionally impose further constraints on the zeros E m of 
i?2«+2 introduced in (3.2) and assume that 

Ex,..., E 2n £ C\{0}, E 2n+ i = 0. (3.4) 

We compactify JC n by adding two points at infinity, Poo + , Poo_, with ^ 
-Poo_, still denoting its projective closure by K, n - Hence K, n becomes a two-sheeted 
Riemann surface of arithmetic genus n. Points P on K. n \{P 00± } are denoted by P — 
(z, y), where y( ■ ) denotes the meromorphic function on tC n satisfying T n {z, y) = 0. 



6 



F. GESZTESY AND H. HOLDEN 



For notational simplicity we will usually tacitly assume that n e N (the case 
n = being trivial) . 

In the following the roots of the polynomials F n and H n will play a special role 
and hence we introduce on C x E 

n n 

F n (z,x) = Y[(z - Hj(x)), H n (z,x) = h (x) ]J(z - Vj{x)), (3.5) 

3=1 3=1 

temporarily assuming 

h o {x)^0, xei (3.6) 

Moreover, we introduce 

p,j(x) = (fij(x), -Hj(x)G„(fi,j(x),x)) e /C„, j = 1, . . . ,n, x £ R, (3.7) 
i>j(a;) = (i/j(a;),i>j (a;), a;)) e /C„, j = 1, . . . , n, x e E, (3.8) 

and 

Po = (0,0). (3.9) 
The branch of y( • ) near Poo ± is fixed according to 

lim , m J ( r? m x = Tl- (3.10) 

|*(P)|-«x> z(P)G n {z(P),x) 

Due to assumption (3.1), u is smooth and bounded, and hence F n (z, ■ ) and H n (z 7 ■ ) 
share the same property. Thus, one concludes 

N ,u k e C(E), j,k = l,...,n, (3.11) 

taking multiplicities (and appropriate reordering) of the zeros of F n and H n into 
account. 

Next, we introduce the fundamental meromorphic function <p( ■ , x) on JC n by 

ai-d \ y-zG n (z,x) zH n (z,x) 

mx)= F n (z,x) = y + zG n {z, X ) > P =^v)^n,xGR. (3.12) 

Assuming (3.4) and (3.6), the divisor (<f>( ■ , x)) of 4>{ ■ , x) is given by 

(<f)(-, x)) =Vp z( x) -V Poo+ ^ x) , (3.13) 

taking into account (3.10). Here we abbreviated 

/£ = {/ti, . . . , /x„}, v = . . . , n } e cr n K n , (3-14) 

where cr m /C„, m e N, denotes the mth symmetric product of lC n . If h is permitted 
to vanish at a point Xi G R, then for a; = Xi, the polynomial 7J„( • , Xi) is at most 
of degree n — 1 (cf. (2.18)) and (3.13) is altered to 

(</>(•, a;i)) =I>p p oo _i> 1 (x 1 ),...,p n _ 1 (x 1 ) ( 3 -15) 

that is, one of the 0j(x) tends to Poc_ as x — > a;i (cf. also (3.36)). Analogously 
one can discuss the case of several z>j approaching Poo_ • Since this can be viewed 
as a limiting case of (3.13), we will henceforth not particularly distinguish the case 
ho ^ from the more general situation where ho is permitted to vanish. 
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Given (j>(- ,x), one defines the associated Baker- Akhiezer vector ^(-,x,Xo) on 
^n\{Poo+ ) Poo^Po} by 

nP^,x )=(^^ x x ^, Pe/C„\{Poo + ,Pco_,Po}, (x,x ) el 2 , (3.16) 
where 

M p ^,x )=e3q?(-{l/z)J dx 1 4>{P,x') - {x - xo^j , (3.17) 

MP,x,x ) = -MP,x,x )<p(P,x)/z. (3.18) 

Although ^ is formally the analog of the Baker-Akhiezer vector of the stationary 
CH hierarchy when compared to analogous definitions in the context of the KdV or 
AKNS hierarchies, its actual properties in a neighborhood of its essential singularity 
feature characteristic differences to standard Baker Akhiezer vectors as discussed 
in [13, Ch. 5] and [12]. 

The basic properties of <j> and 'J then read as follows. 

Lemma 3.1. Suppose (3.1), assume the nth stationary CH equation (2.22) holds, 
and let P = (z,y) e IC n \{P 00+ ,P oc _,Po}, (x,Xo) <E K 2 - Then <fi satisfies the 
Riccati-type equation 



MP, x) - z-^iP, xf - 2<j>{P, x) + Au(x) - u xx (x) = 0, 

as well as 

zH n (z, x) 



<P{P,x)cj>{P\x) 
4>{P,x)+4>{P*,x)= -2 

4>{P,x)-4>(P*,x) 



F n (z,x) ' 
zG n (z, x) 
F n (z,x) ' 
2y 



F n {z,xY 
while ^ fulfills 

V x (P,x,x ) = U{z,x)V{P,x,x ), 

- y*(P, x, x ) = zV n (z, ar)*(P, x, x ), 

MP,x,xo)=( ~i Z,X \ ) 1 expf-iy/z) f dx'F n (z, x')' 1 
\F n {z,x )J \ J X0 

F n {z,x) 



ip 1 (P,x,x )ipi(P*,x,x ) = 

1p2(P,X,X )tp 2 (P*,X,X ) = 



F n (z,x ) ' 
H n (z,x) 
zF n (z,x ) ' 

G n (z, x) 



ip 1 (P,x,x )i)2{P*,x,x Q ) + tpi(P*,x,x a )^ 2 (P, x,x ) = 2 
ip 1 (P,x,x )ip 2 (P*,x,x ) - ipi{P* , x, x )il> 2 (P, x, x ) 



F n (z,x ) ' 
22/ 



zF n (z,x ) ' 



3.19) 

3.20) 
3.21) 
3.22) 

3.23) 
3.24) 

3.25) 

3.26) 

3.27) 

3.28) 

3.29) 



In addition, as long as the zeros of F n ( ■ , x) are all simple for I £ SI, fl C I an 
open interval, <!'(• , x,xq), x,x n e ft, is meromorphic on /C„\{Pq}. 
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Next, we recall the Dubrovin-type equations for Since in the remainder of 
this section we will frequently assume K n to be nonsingular, we list all restrictions 
on IC n in this case, 

E m G C\{0}, E m ^ E m , for m j= m', m, m' = 0, . . . , 2n + 1, E 2n+1 = 0. (3.30) 

Lemma 3.2. Suppose (3.1) and the nth stationary CH equation (2.22) holds subject 
to the constraint (3.30) on an open interval fl^ C M. Moreover, suppose that 
the zeros [ij, j = l,...,n, of F n (-) remain distinct and nonzero on f2 M . Then 
{fij}j=i,..., n , defined by (3.7), satisfies the following first-order system of differential 
equations 

fi jlX (x) = 2 V ^j X } ] flfaix) - »e(x))-\ j = l,...,n,xG%. (3.31) 

Next, assume }C n to be nonsingular and introduce the initial condition 

{fij(xo)}j=i,...,n C K n (3.32) 

for some x 6 K, where Hj(xo) ^ 0, j = 1, . . . , n, are assumed to be distinct. Then 
there exists an open interval fi M C M, with x n E 0^, such that the initial value 
problem (3.31), (3.32) has a unique solution {fij}j=i,... in C fC n satisfying 

A, eC°°(0„U i = l,...,n, (3.33) 

and fjij, j = 1, . . . ,n, remain distinct and nonzero on Cl^. 

Combining the polynomial approach in Section 2 with (3.5) yields trace formulas 
for the CH invariants. For simplicity we just record two simple cases. 

Lemma 3.3. Suppose (3.1), assume the nth stationary CH equation (2.22) holds, 
and let Then 

_^ n _^ 2n+l 

j—l rn — 

4u(a0-u xx (a0 = -( fl E ™ j ( II ^)~ 2 J • (3-35) 

\ m=0 / \j=l / 

Next we turn to asymptotic properties of (j) and V'j, j = 1, 2. 

Lemma 3.4. Suppose (3.1), assume the nth stationary CH equation (2.22) holds, 
and let P = (z,y) G /C„\{Poo + , P 00 _ , P }, z G M. T/ien 

(p , = f-2C- 1 -2«(*)+ «»(*)+ O(C), _r 
rV ' 'c-o|2u(a;) + «*(*)+ 0(C), P^Poc_, ' 

(2n+l \ V 2 

J] /nW-'C + O^ 2 ), P^Po, C = z 1/2 , (3.37) 
m=0 / 

and 

MP,x,x ) = o exp(±(a:-a;o))(l + 0(0)» P ^ P°o±, C = V«, (3-38) 
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MP , X , X0) c = exp( ± (,-,o))| (2 2 u ;^ W)c+0(c ^ p p ZpZ (3 - 39) 

C = l/z, 

1/2 



MP,x,x ) ^cxp^-i jf rfx'^ [] ^^O-'+Oa)^, (3.40) 



P - P , C = ^ 1/2 , 

^ 2 (P, a; , a; o) c = o O(r 1 )exp^-i^dx'^ jf] MxT 1 + 0(1)^, 

P^Po, C^z 1 ' 2 . (3.41) 

Since the representations of <fr and u in terms of the Riemann theta function 
associated with K n (assuming K n to be nonsingular) are not explicitly needed in 
this paper, we omit the corresponding details and refer to the detailed treatment 
in [13, Ch. 5] and [12] instead. 

Finally, we will recall that solvability of the Dubrovin equations (3.31) onD^Ci 
in fact implies equation (2.22) on 

Theorem 3.5. Fix n e N, assume (3.30), and suppose that {{ij}j=i,...,n satisfies 
the stationary Dubrovin equations (3.31) on an open interval f2 M C R such that 
j = 1, ... ,n, remain distinct and nonzero on Cl^. Then u e defined by 

< x ) = 2 E w w-^E ^ ( 3 - 42 ) 

j= 1 m=0 

satisfies the nth stationary CH equation (2.22), that is, 

s-CH„(w) = on %. (3.43) 

4. Basic Facts on Self-adjoint Hamiltonian Systems 

We now turn to the Weyl-Titchmarsh theory for singular Hamiltonian (canon- 
ical) systems and briefly recall the basic material needed in the following section. 
This material is standard and can be found, for instance, in [6], [16], [17], [18], [19], 
[20], and the references therein. 

Hypothesis 4.1. (i) Define the 2x2 matrix J — (" and suppose aj.k,bj_k 6 
AocW> 3,k = 1,2 and A(x) = (a jtk (x)) . fe=12 > 0, B(x) = (&,,fc(^)^ =1 , 2 = 
B{x)* for a.e. x G M. We consider the Hamiltonian system 

J&'(z,x) = {zA(x)+B(x))&(z,x), zeC (4.1) 

for a.e. x e M, where z plays the role of the spectral parameter, and where 

&(z,x) = (M*,x) Mz,x)) T , <Pj(z, •) e ACloc(K), j = 1,2. (4.2) 

//ere j4Ci c(R) denotes the set of locally absolutely continuous functions on R arte? 
t/ie M* and M T denote the adjoint and transpose of a matrix M, respectively. 
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(ii) For all nontrivial solutions \P of (4.1) we assume the deftniteness hypothesis 
(cf. [4, Sect. 9.1]) 

r d 

dx &(z, x)*A(x)&(z, x)>0, (4.3) 



/ 

J c 



on every interval (c, d) C R, c < d. 

A simple example of a Hamiltonian system satisfying (4.3) is obtained when 
A(x) = ( w ^ ° ) , for some weight function w G i/^ )c (]R), w > a.e. on R, and 
^2,2(2) > a.e. on R (cf. Section 5). Hypothesis 4.1 (ii) clearly holds in this case. 

Next, we introduce the vector space (—00 < a < b < 00) 

L A ((a,b)) = |(/> : (a, b) — ► C 2 measurable J dx (<fi(x), A(x)(j)(x)) C 2 < oo|, 

(4.4) 

where ((f), V)c 2 = Sj=i ^jV'j denotes the standard scalar product in C 2 . Fix a point 
xq G R. Then the Hamiltonian system (4.1) is said to be in the limit point case at 00 
(resp., —00) if for some (and hence for all) z € C\K, precisely one solution of (4.1) 
lies in L 2 A ((x , 00)) (resp., L A ((—oo,x ))). (By the analog of Weyl's alternative, if 
(4.1) is not in the limit point case at ±00, all solutions of (4.1) lie in L A ((x , ±00)) 
for all z G C. In the latter case the Hamiltonian system (4.1) is said to be in the 
limit circle case at ±00.) 

To simplify matters for the remainder of this section, we will always suppose the 
limit point case at ±00 from now on. 

Hypothesis 4.2. Assume Hypothesis 4-1 and suppose that the Hamiltonian system 
(4.1) is in the limit point case at ±00. 

An elementary example of a Hamiltonian system satisfying Hypothesis 4.2 is 
given by the case where all entries of A and B are essentially bounded on R (cf. 
Section 5). 

When considering the Hamiltonian system (4.1) on the half-line [20,00) (resp., 
(— 00, xo]), a self-adjoint (separated) boundary condition at the point x is of the 
type 

a^(x ) = 0, (4.5) 

where a = (ol\ a.2) G C lx2 satisfies 

eta* = I, aJa* = (equivalently, |ai| 2 + |«2| 2 = 1, Im(a2oT) = 0). (4-6) 

In particular, the boundary condition (4.5) (with a satisfying (4.6)) is equivalent 
to aiipi(xo) + ct^V^^o) = with ai/a 2 G R if a 2 ^ and a^/ai G R if ct\ ^ 0. 
The special case cto = (1 0) will be of particular relevance in Section 5. Due to our 
limit point assumption at ±00 in Hypothesis 4.2, no additional boundary condition 
at ±00 needs to be introduced when considering (4.1) on the half-lines [x ,oo) 
and (— 00, a;o]. The resulting full-line and half- line Hamiltonian systems are said 
to be self-adjoint on R, [x , 00), and (— 00, £o] ; respectively (assuming of course a 
boundary condition of the type (4.5) in the two half- line cases). 

Next we digress a bit and briefly turn to Hcrglotz functions and their representa- 
tions in terms of measures, the focal point of Wcyl Titchmarsh theory (and hence 
spectral theory) of self-adjoint Hamiltonian systems. 
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Definition 4.3. Any analytic map m: C+ — > C+ is called a Herglotz function (here 
C+ = {z G C | Im(z) > 0}). Similarly, any analytic map M : C+ -> C fcxfc , k G N, is 
called a k x k matrix-valued Herglotz function if Im(M(z)) > for all z G C+. 

Herglotz functions are characterized by a representation of the form 

/oo 
duj(\)((\- z)- 1 - Xil + X 2 )- 1 ), zgC\R, (4.7) 
-OO 

/OO 
duj(X) (1 + A 2 ) -1 < oo, (4.8) 
-oo 

w((A, n]) = limlim — / aV Im (m(z/ + ie)) , (4.9) 

5i0 ej.0 7T J x+S 

in the following sense: Every Herglotz function admits a representation of the 
type (4.7), (4.8) and conversely, any function of the type (4.7), (4.8) is a Herglotz 
function. Moreover, local singularities and zeros of m are necessarily located on the 
real axis and at most of first order in the sense that 

= lim(w(A + e) - w(A - e)) = - \imie m(X + ie) > 0, A G R, (4.10) 

ej.0 ej.0 

limiemiX + iey 1 > 0, A G R. (4.11) 

ej.0 

In particular, isolated poles of m are simple and located on the real axis, the 
corresponding residues being negative. Analogous results hold for matrix-valued 
Herglotz functions (see, e.g., [14] and the literature cited therein). 

For subsequent purpose in Section 5 we also note that — 1/z is a Herglotz function 
and compositions of Herglotz functions remain Herglotz functions. In addition, 
diagonal elements of a matrix-valued Herglotz function are Herglotz functions. 

Returning to Haniiltonian systems on the half- lines [xq, ±oo) satisfying Hypothe- 
ses 4.1 and 4.2, we now denote by ^±(z, x, xq) the unique solution of (4.1) satisfying 
i $>±(z, - ,xo) G L\{[xq, ±oo)), z G C\M, normalized by ipi,±(z,xo,xo) = 1. Then 
the half-line Weyl-Titchmarsh function m±(z,x), associated with the Hamiltonian 
system (4.1) on [x, ±oo) and the fixed boundary condition ao = (1 0) at the point 
x G R, is defined by 

m±(z, x) = ip2,±(z,x, xo)/ipi, ±(z,x,xo), z G C\R, ±x > xq. (4-12) 

The actual normalization of $f±(z,x,Xo) was chosen for convenience only and is 
clearly immaterial in the definition of m±(z, x) in (4.12). 

One easily verifies that m±(z,x) satisfies the following Riccati-type differential 
equation, 

m'(z, x) + [62,2(3;) + a,2,2(x)z]m(z, x) 2 (4-13) 
+ [61,2(2;) + 62,1(2;) + (01,2(0;) + d2,i(x))z]m(z,x) + 61,1(2;) + ai t i(x)z = 0. 

Finally, the 2x2 Weyl-Titchmarsh matrix M(z,x) associated with the Hamil- 
tonian system (4.1) on R is then defined in terms of the half-line Weyl-Titchmarsh 
functions m±(z,x) by 

M(z, x) = (Mjj, (z, x)). Jl=lt2 , z G C\R, (4.14) 

Mi.i(z, x) — [m,-(z, x) — m+(z, a;)] -1 , 

M li2 (z,x) = M 2 ,i(z,x) 
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= 2~ 1 {m-(z,x) - m + (z, x)]~ 1 [m-(z, x) + m+(z,x)], (4-15) 

M 2 ^ 2 (z,x) = [m-(z,x) — m + (z,x)]~ 1 m-(z,x)m + (z,x). 

One verifies that M(z,x) is a 2 x 2 matrix- valued Herglotz function. We emphasize 
that for any fixed xq G R, M(z,xq) contains all the spectral information of the 
sclf-adjoint Hamiltonian system (4.1) on R (assuming Hypotheses 4.1 and 4.2). 

5. Real- Valued Algebro-geometric CH Solutions 
and the Associated Isospectral Torus 

In our final and principal section we study real-valued algebro-geometric solu- 
tions of the CH hierarchy associated with curves K, n whose affine part is nonsingular 
and determine the isospectral manifold of smooth bounded CH solutions. We focus 
on the stationary case as this is the primary concern in this context. 

To study the direct spectral problem we first introduce the following assumptions. 

Hypothesis 5.1. Suppose 

E <E 1 <---<E 2n < E 2n+1 = (5.1) 
and let u be a real-valued solution of the nth stationary CH equation (2.22), 

s-CH„(u) = 0, (5.2) 
(i.e., u is a particular algebro-geometric CH potential), satisfying 

ueC°°(R), d^ue L°°(K.), fc = 0,1,2, (5.3) 
4u - u xx > 0. (5.4) 

We start by noticing that the basic stationary equation (3.23), 

* x (z, x) = U(z, x)V{z, x), * = (V>i, ^2) T , (z, i)eCxI. (5.5) 
is equivalent to the following Hamiltonian (canonical) system 

J* x (z,x) = [zA(x)+B(x)]*(z,x), $ = (^i,^ 2 ) T , (z,x) e C x R, (5.6) 
where 

J=(° 1 "q 1 ), *{z,x) = *(z,x), z = -1/z, (5.7) 

A{x) = ^(x)^u xx{ x) B{x) = ^ -iy x£R _ (58) 

In particular, due to assumptions (5.3) and (5.4), the Hamiltonian system (5.6) 
satisfies Hypotheses 4.1 and 4.2. Explicitly, the Hamiltonian system (5.6) boils 
down to 

i> 1>x (z,x)=i} 2 (z,x)-i}i(z,x), (5.9) 

4>2,x(z, x) = -z(4u(x) - u xx (x))^!(z,x) + 4> 2 (z,x), (z,i)eCxI, (5.10) 

and upon eliminating ?/>2 results in a particular case of the weighted Sturm-Liouville 
problem 

1 

r 

of the type 

-i>i,xx{z,x) + $i(z,x) = z(4u(x) - u xx (x))i}i(z,x), (z,x) G C x R, (5.12) 



dx^ dx ^ 



(5.11) 
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X < 



A e M, (5.14) 



with "weight" r = (4m — u xx ) and constant coefficients p = q = 1. 
Introducing 

n 

X={J[E 2l ,E 2e+1 ], (5.13) 

we define 

i? 2 „+ 2 (A) 1 / 2 = |i? 2 „ +2 (A) 1 / 2 | 

-1 for A e (S 2n+ i,oo), 

(-1)" +J for A e (E 2j _ 1 ,E 2j ), j = l,...,n, 
(-1)" for A e (— oo, S ), 

for A g {E2jjE2j+i)j j = 0, . . . , n, 

and 

i?2n+ 2 (A) 1/2 = limi? 2 „ +2 (A + 2e) 1 /2 ; AGS, (5.15) 

and analytically continue R 2n +2{ • )^ 2 to C\S. We also note the property 

iWp=4 +2 W 1/2 (5.16) 

For notational convenience we will occasionally call (E 2 j_\,E 2 j), j = l,...,n, 
spectral gaps and E 2 j_\,E 2 j the corresponding spectral gap endpoints. 
Next, we introduce the cut plane 

n = C\E, (5.17) 
and the upper, respectively, lower sheets II-i- of K, n by 

n± = {(z, ±R 2n+2 {zf' 2 ) eK, n \zeU} (5.18) 

with the associated charts 

C±:II ± ^n, P={z,±R 2n+2 (z) 1 ' 2 )^ z. (5.19) 

The two branches ^±(z,x,x ) of the Baker-Akhiezer vector ^(P,x,x ) in (3.16) 
are then given by 

* ± (z,x,xo) = *(P,x,x ), P=(z,y)eU ± , tf± = (V>i,±, V>2,±) T , (5.20) 
and one infers from (3.38) that 

i>i,±(z, ■ ,xo) £ L 2 ((x , Too)) for \z\ sufficiently large. (5.21) 
Thus, introducing 

*±(z,x,x ) = V T (z,x,x ), *± = (Va,±,V>2,±) T , z = -l/z, (5.22) 

and the two branches </>±(z, x) of </>(P, x) on II ± by 

<t>±{z,x) = <t>(P,x), P={z,y)eU ± , (5.23) 

one infers from (4.12) and (5.21) that the Weyl-Titchmarsh functions m±(z,x) 
associated with the self-adjoint Hamiltonian system (5.6) on the half-lines [x, ±oo) 
and the Dirichlet boundary condition indexed by ceo = (1 0) at the point x G M, 
are given by 

rh±(z,x) = 4>2,±{z,x,x )/ip li ±(z,x,xo) = ip 2 , T (z,x,x )/ipi, T {z,x,x ) 

= (-l/z)<j> T (z,x), zgC\S. (5.24) 



14 



F. GESZTESY AND H. HOLDEN 



More precisely, (5.21) yields (5.24) only for \z\ sufficiently large. However, since 
by general principles m±(-,x) are analytic in C\K, and by (3.12), <fi±(-,x) are 
analytic in C\S, one infers (5.24) by analytic continuation. In particular, (5.21) 
extends to all z E C\S, that is, 

V>i,±(z> -,xo) G L 2 ((a;o,T^)), z G C\S. (5.25) 

Next, we mention a useful fact concerning a special class of Hcrglotz functions 
closely related to the problem at hand. The result must be well-known to experts, 
but since we could not quickly locate a proof in the literature, we provide the simple 
contour integration argument below. 

1 12 

Lemma 5.2. Let Pn be a monic polynomial of degree N. Then PN/R^n+z * s a 
Herglotz function if and only if one of the following alternatives applies: 
(i) N = n and 

n 

Pn(z) = Y[{z- aj ), aj E [E 2j -i,E 2j ], j = l,...,n. (5.26) 
If (5.26) is satisfied, then P n /R\n+2 admits the Herglotz representation 

MX) 1 f \Pn(X)\dX 1 

Rm + 2{zy/*- Kh\R2 n+ i{XmX-z' ZG(L ^- (5 ' 27j 
(u) N = n-\-l and 

n 

P n +i{z) = X\{z-b i ), bo € (—oo, Eq], bj G [E 2j -i,E 2j ], j = 1, ...,n. (5.28) 

£=0 

1 /2 

// (5.28) is satisfied, then P n+ i/R 2 ' n+2 admits the Herglotz representation 

Pn+i(z) = ( P n+ l(l) \ If \P n+ i(X)\dX ( 1 A 

i?2n +2 (z)V2 U2n+2W 1/2 y 7T 7 E |i? 2n+2 (A)l/2 | ^ A - Z 1 + A 2 , 

z G C\E. (5.29) 

Proof. Since Herglotz functions are O(z) as |z| — > oo and cannot vanish faster than 
0(1/ z) as |z| — > oo, we can confine ourselves to the range N G {n, n+ 1, n + 2}. We 
start with the case N — n and employ the following contour integration approach. 
Consider a closed oriented contour r#. e which consists of the clockwise oriented 
semicircle C £ — {z e C | z = Eo—ecxp(—ia), —n/2 < a < ir/2} centered at Eq, the 
straight line L + = {z E C+ | z = Eo+x+ie, < x < R} (oriented from left to right), 
the following part of the counterclockwise oriented circle of radius (R 2 + e 2 ) 1 / 2 
centered at E , C R = {z e C | z = E + (R 2 + e 2 ) 1 / 2 exp(i/3), arctan(e/i?) < (3 < 
2ir — arctan(e/i?)}, and the straight line L_ = {z E C_ | z = E +x — ie, < x < R} 
(oriented from right to left). Then, for e > small enough and R > sufficiently 
large, one infers 

Pn(z) 1 / 1 P„(C) 



R2n+2(zy/ 2 2m Jr R E C - z i? 2 „+ 2 (C) 1 / 2 

= 1 f 1 PnWdX 

XRtoo tt y E A - z iR 2n+2 (Xy/ 2 



sio, 
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Here we used (5.14) to compute the contributions of the contour integral along 
[E , R] in the limit e [ and note that the integral over Cr tends to zero as R | oo 
since 

idlW^ 0<lcrl) - (5 - 31) 

Next, utilizing the fact that P n is monic and using (5.14) again, one infers that 
F n (X)dX/[iR 2n+2 (X) 1 ^ 2 ] represents a positive measure supported on £ if and only 
if P n has precisely one zero in each of the intervals {E 2 j~i, E 2 j], j — 1, . . . ,n. In 
other words, 

if and only if P n has precisely one zero in each of the intervals [E2j-i, E 2 j], j = 
1, . . . ,n. The Hcrglotz representation (4.7), (4.8) then finishes the proof of (5.27). 

In the case where N = n+1, the proof of (5.28) follows along similar lines taking 
into account the additional residues at ±i inside Tr^ which are responsible for the 
real part on the right-hand side of (5.29). 

Finally, in the case N = n + 2, assume that P n +2/R 2 n+2 i s a Herglotz function. 
Then necessarily, 

D Pn+2 }"] /2 =a + bz+ f° dco(X) (A - z)-\ z G C\E (5.33) 

R2n+2(Z) 1 ' 2 Jeo 

for some a G M, b > 0, and some finite (positive) measure u> supported on [i?o,0], 
since 

limIm(P„ +2 (A)i? 2 „ +2 (A + is)- 1/2 ) = for A > E 2n+2 = and A < E . (5.34) 

In particular, (5.33) implies 

Pn+2(z)R 2n +2(z)- 1/2 : ^ bz + 0(l), b>0. (5.35) 

However, by (5.14), one immediately infers 

j P„+ 2 (A)i 1 , 2n+2(A)- 1 / 2 = -A + 0(1). (5.36) 

Afoo 

This contradiction dispenses with the case N = n + 2. □ 

Now we are in position to state the following result concerning the half-line 
and full-line Weyl-Titchmarsh functions associated with the self-adjoint Hamil- 
tonian system (5.6). We denote by fh±(z, x) the Weyl-Titchmarsh m-functions 
corresponding to (5.6) associated with the half-lines (x, ±oo) and the Dirichlet 
boundary condition indexed by ao = (1 0) at the point i£l, and by M(z,x) the 
2x2 Weyl-Titchmarsh matrix corresponding to (5.6) on K (cf. (4.12), (4.14), and 
(4.15)). Moreover, S° denotes the open interior of S and the real part of a matrix 
M is defined as usual by Re(M) = (M + M*)/2. 

Theorem 5.3. Assume Hypothesis 5.1 and let (z,x) G M x (C\S), z = —1/z. 
Then 
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1±R / ^2n+ 2 (i) 1/2 \ ^G n (jij(x),x)(lTej(x)) 



\ iF n {i,x) J dF n ([i :j (x),x)/dz z - fij(x) 

, 1 f \R2n + 2{\) 1/2 \d\ ( 1 A_\ 

± 7 ri s |AF„(A,x)| \\-z \ + 
where £j(x) e {1, —1}, i = 1, • • • ,n, is chosen such that 
G n (nj(x),x)£j(x) 



dF n (/j,j(x),x)/dz 

Moreover, 



>0, j = l,...,n. (5.39) 



M(~ z , x) = =L_ Hf»<*' s > *^<*' f\ (5.40) 

1 ; 2i? 2n+2 (z) 1 /2 ^zG„(z,a;) zF n (z,x) J v y 

= Re(M(i, *)) + £ dfi(A, x) (^i- - , (5.41) 



where 



0(A,,) = . p 1 ml/a f -f" ( ^|l A e E°. (5.42) 

27r«i? 2 „+2(A) 1 /2 \ v -AG„(A,x) -AF„(A,x)/' 

TTie essential spectrum of the half-line Hamiltonian systems (5.6) on [x, ±oo) (i/rai/i 
any self-adjoint boundary condition at x) as well as the essential spectrum of the 
Hamiltonian system (5.6) on K is purely absolutely continuous and given by 

n-1 

|J [ - E%\ -E£ +1 ] U [ - E£ , oo) . (5.43) 

e=o 

The spectral multiplicities are simple in the half-line cases and of uniform multi- 
plicity two in the full-line case. 

Proof. Equation (5.37) follows from (3.12), (5.14), and (5.24). Equation (5.40) is 
then a consequence of (3.20)^(3.22), (4.14), (4.15), (5.24), and (5.37). Different 
self-adjoint boundary conditions at the point x lead to different half-line Hamil- 
tonian systems whose Weyl-Titchmarsh functions are related by a linear fractional 
transformation (cf., e.g., [6]), which leads to the invariance of the essential spec- 
trum with respect to the boundary condition at x. In order to prove the Herglotz 
representation (5.38) one can follow the corresponding computation for Schrodingcr 
operators with algebro-geometric potentials in [21, Sect. 8.1]. For this purpose one 
first notes that by (5.29) also R2 n +2{z) 1 ^ 2 /[zF n {z, x)] is a Herglotz function. A 
contour integration as in the proof of Lemma 5.2 then proves 

iW2(*) 1/2 _ R j R2n + 2(i) 1/2 \ , J- l^n+2( W W) 1/2 | 1 



zF n (z,x) \ iF n (i,x) J Hj(x)\dF n (fj,j(x),x)/dz\ z - fij(x) 

, 1 f \R2n+2W 1/2 \d\f 1 A 



tt7 s |af„(a,z)| Va-^ i + a 2 

Re / R2n+2(i) 1/2 \ _ y- G n (Hj(x),x)Sj{x) 

V iF n (i,x) ) 4^ dF n (fj,j(x),x)/dz 



(5.44) 



V iF n (i,x) ) dF n {nj{x),x)/dz z- Hj{x) 

i r \R 2n+2 (\md\ ( i a_\ 

+ W S \XF n (X,x)\ {X-z 1 + A 2 / °> 
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The only difference compared to the corresponding argument in the proof of Lemma 
5.2 concerns additional (approximate) semicircles of radius e centered at each Hj(x), 
j = 1, ...,n, in the upper and lower complex half-planes. Whenever, fj,j(x) G 
(E 2 j-i, E 2 j), the limit £ J. picks up a residue contribution displayed in the sum 
over j in (5.44). This contribution vanishes, however, if fJ,j(x) € {E2j-i, E%j}. In 
this case dF n (fij(x), x)/dz ^ by (4.10) and R.2 n +2(tij(x)) = by (2.17). Equation 
(5.45) then follows from (3.7) and the sign of sj(x) must be chosen according to 
(5.39) in order to guarantee nonpositive residues in (5.45) (cf. (4.10)). 

Next, we apply the Lagrange interpolation formula. If Q n -i is a polynomial of 
degree n — 1, then 

Qn-i(z) = F n ( Z )j2 , Q P n ; l{ ^] — — ' zeC - ( 5 - 46 ) 

dF n {m)ldzz-Hi 

Since F n and G n are monic polynomials of degree n, we can apply (5.46) to Q n -i = 
G n — F n and obtain 

Gn(z,x) =1 | y G n (^(x),x) 1 

F n (z,x) ^dF n (fj,j(x),x)/dz z- fj,j{x)' 

Insertion of (5.47) into (5.37) then yields 

rh±(z, x) = ±R ^( Z ) 1/2 + y G n ( j i J (x),x)[s 1 (x) + (l-e 1 (x))] 1 + 
zF n (z,x) j-^ dF n (fij(x),x)/dz z - /j,j(x) 

(5.48) 

and hence (5.38) follows by inserting (5.45) into (5.48). Equations (5.41) and (5.42) 
are clear from the matrix analog of (4.9). 

The statement (5.43) for the essential half-line spectra then follows from the fact 
that the measure in the Herglotz representation (5.38) of m± (as a function of z) 
is supported on the set £ in (5.13), with a strictly positive density on the open 
interior E° of E. The transformation z — > —1/z then yields (5.43) and since half- 
line spectra with a regular endpoint x have always simple spectra this completes 
the proof of our half-line spectral assertions. The full-line case follows in exactly 
the same manner since the corresponding 2x2 matrix-valued measure f2 in the 
Herglotz representation (5.41) of M (as a function of z) also has support S and 
rank equal to two on S°. □ 

Returning to direct spectral theory, we note that the two spectral problems (5.6) 
on M associated with the vanishing of the first and second component of ^ at x, 
respectively, are clearly self-adjoint since they correspond to the choices a = (1 0) 
and a = (0 1) in (4.5). Hence a comparison with (3.5), (3.26), and (3.27) necessarily 
yields {tij(x)}j=i t ... in , { u j( x )}j=i,---,n C K. Thus we will assume the convenient 
eigenvalue orderings 

Hj(x) < /j,j + i(x), Vj{x) < Vj + i(x) for j = 1, ... ,n — 1, x£l. (5.49) 

The zeros of ipi ( • , x) belong to the Dirichlet spectral problem associated with the 
Hamiltonian system (5.6) (resp., the weighted Sturm-Liouville problem (5.12)) on 
R. A comparison with (3.26) then relates the zeros fij(xi), j = 1, . . . , n, of F n ( ■ , x\) 
in (3.5) to the Dirichlet spectrum of (5.6) (resp., (5.12)) on R. The correspondence 
between each /x,- and the related spectral point of the Dirichlet problem (5.6) (resp., 



18 



F. GESZTESY AND H. HOLDEN 



(5.12)) on R is of course effected by the transformation z — > —1/z. In contrast to 
this, the zeros of 4> 2 (- , Xi) do not belong to the Neumann spectrum associated 
with the Hamiltonian system (5.6) (resp., the weighted Sturm-Liouville problem 
(5.12)) on R. In fact, by (5.9), zeros of ip 2 ( ■ , x\) correspond to a mixed boundary 
condition of the type ifi, x (xi) + ipi(xi) = 0. By (3.27), this relates the zeros 
Vj(x\), j = 1, ...,n, of H n (- ,x\) in (3.5) to the spectrum of (5.6) (resp., (5.12)) 
on R corresponding to the self-adjoint boundary condition ip i tX {%i) + ipi(xi) = 0. 

Combining Lemma 5.2 with the Herglotz property of the 2x2 Weyl-Titchmarsh 
matrix M( ■ , x) then yields the following refinement of Theorem 3.2. 

Theorem 5.4. Assume Hypothesis 5.1. Then {Aj}j=i,...,n> with the projections 
Hj(x), j = 1, . . . ,n, the zeros of F n (- ,x) in (3.5), satisfies the first-order system of 
differential equations (3.31) on O p = R and 

fij eC°°(R,/C„), j = l,...,n. (5.50) 

Moreover, 

Hj{x) G [E 2j -i,E 2j ], j = 1, . . . , n, x G R. (5.51) 

In particular, fij(x) changes sheets whenever it hits E 2 j-i or E 2 j and its projection 
Hj(x) remains trapped in [£27-1, E 2 j] for all j = 1, . . . , n and x G R. The analogous 
statements apply to Vj(x) and one infers 

Vj{x) e [E 2 ^ u E 2j ], j = 1, . . . ,71, x e R. (5.52) 

Proof. Since M( • , x) is a 2 x 2 Herglotz matrix, its diagonal elements are Herglotz 
functions. Thus, 

M hl (z,x)= Hn{z lf M^z,x)= ~ XFn{ *l*L , S=-l/z (5-53) 
2R 2n+2 (z) L / z 2R 2n+2 (z) L / z 

are Herglotz functions (the left-hand sides with respect to z, the right-hand sides 
with respect to z) and the interlacing properties (5.51), (5.52) then follow from 
(5.28) and (5.26). □ 

Remark 5.5. Combining the interlacing property (5.51) with (2.18), (2.19), and 
(2.20) yields (cf. also (3.35)) 

Au(x) -u xx {x) = -( ]J ^m) (n^W 2 ) > 0, xeR, (5.54) 

^m=0 ' ^ j=l ' 

in accordance with (5.4). Moreover, since by (5.52) the Vj(x) also remain trapped 
in the intervals [E 2 j_i, E 2 f\ for all none of the Vj can reach P oc _ and hence 

h = 4m + 2u x ^= on R (cf. the discussion surrounding (3.15)). Actually, 

ho(x) > 0, lei, (5.55) 

1/2 

since H n ( ■ , x)/R 2n+2 is a Herglotz function (cf. (5.27)). 

Remark 5.6. The zeros fij(x) G (E 2 j_i, E 2 j), j = 1, . . . , n of F n ( ■ , x) which are 
related to eigenvalues of the Hamiltonian system (5.6) on R associated with the 
boundary condition tpi(x) — 0, in fact, are related to left and right half-line eigenval- 
ues of the corresponding Hamiltonian system restricted to the half-lines (— 00, x] and 
[x, 00), respectively. Indeed, by (5.22) and (5.25), depending on whether fij(x) G n+ 
or fij(x) G n_, jUj(x) is related to a left or right half-line eigenvalue associated with 
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the Dirichlet boundary condition r/>i(ar) = 0. A careful investigation of the sign of 
the right-hand sides of the Dubrovin equations (3.30) (combining (5.1), (5.14), and 
(5.18)), then proves that the jJbj(x) related to right (resp., left) half-line eigenvalues 
of the Hamiltonian system (5.6) associated with the boundary condition tpi (x) = 0, 
are strictly monotone increasing (resp., decreasing) with respect to x, as long as 
the Hj stay away from the right (resp., left) endpoints of the corresponding spec- 
tral gaps (E 2 j-i,E2j). Here we purposely avoided the limiting case where some of 
the Hk(x) hit the boundary of the spectral gaps, /ik(x) G {E 2 k-i, E>2k}> since the 
half-line eigenvalue interpretation is lost as there is no L 2 ((x, ±oo)) 2 eigenfunction 
^(x) satisfying tpi(x) — in this case. In fact, whenever an eigenvalue /J,k(x) hits 
a spectral gap endpoint, the associated point fij (x) on K n crosses over from one 
sheet to the other (equivalently, the corresponding left half-line eigenvalue becomes 
a right half-line eigenvalue and vice versa) and accordingly, strictly increasing half- 
line eigenvalues become strictly decreasing half-line eigenvalues and vice versa. In 
particular, using the appropriate local coordinate (z — E 2 k) 1 ^ 2 (resp., (z — E 2 k-i) 1 ^ 2 ) 
near E 2 k (resp., E 2 k-i), one verifies that Hk(x) does not pause at the endpoints E 2 k 
and E 2 k-\. 

Next, wc turn to the inverse spectral problem and determine the isospectral 
manifold of real-valued, smooth, and bounded CH solutions. 
Our basic assumptions then will be the following: 

Hypothesis 5.7. Suppose 

E <E 1 <---<E 2n < E 2n+1 = 0, (5.56) 

fix Xq G R, and assume that the initial data 

{fij(xo) = (p 3 (x ),-tJ.j(x )G n (iJ, j (xo),xo))} j=lt ,.^ n c K n (5.57) 

for the Dubrovin equations (3.31) are constrained by 

fij(x ) S [E 2 j-i,E 2j ], j = l,...,n. (5.58) 

Theorem 5.8. Assume Hypothesis 5.7. Then the Dubrovin initial value problem 
(3.31), (5.57), (5.58) has a unique solution {Aj}j=i,...,n C JC n satisfying 

^ eC°°(R,)C n ), j = l,...,n, (5.59) 

and the projections fij remain trapped in the intervals {E 2 j_i, E 2 j], j = l,...,n, 
for all x G M, 

pL (x) e [E 2j _ u E 2j ], j = 1, . . . ,n, x G R. (5.60) 
Moreover, u defined by the trace formula (3.34), that is, 

^ n 2n+l 

«(a:) = -^ ft ( a ;)--^B m , x e R, (5.61) 

j— 1 m=0 

satisfies Hypothesis 5.1, that is, 

ueC oc (R), u is real-valued, (5.62) 

d k x u g L°°(R), k G N , (5.63) 

Au -u xx > Q, x G R, (5.64) 
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and the nth stationary CH equation 

s-CH„(u) = on R, (5.65) 

with integration constants C£ in (5.65) given by eg = ci(E_), I = 1, . . . ,n, according 
to (2.26), (2.27). 

Proof. Given initial data constrained by /ij(xo) G (E 2 j-\, E>2j), j = 1,. one 

1 /2 

concludes from the Dubrovin equations (3.31) and the sign properties of R 2n+2 011 
the intervals [E 2 k-i, £-2fe] 7 A; = 1, . . . , n, described in (5.14), that the solution Hj{x) 
remains in the interval \E 2 j_\, E 2 j] as long as fij(x) stays away from the branch 
points (E2j-i,0), (E2j,0). In case fij hits such a branch point, one can use the 
local chart around (E m ,0), with local coordinate ( — a(z — Em) 1 / 2 , a G {1, — 1}, 
to e {2j - l,2j}, to verify (5.59) and (5.60). Relations (5.61)-(5.63) are then 
evident from (5.59), (5.60), and 

\d^j{x)\<C k , t6N ,j = l n,i6R. (5.66) 

In the course of the proof of Theorem 3.5 presented in [12] (cf. also [13, Sect. 
5.3]), one constructs the polynomial formalism (F n , G n , H n , R2n+2, etc.) and 
then obtains identity (3.35) as an elementary consequence. The latter immediately 
proves (5.64). Finally, (5.65) follows from Theorem 3.5 (with f2„ = M.). □ 

Corollary 5.9. Fix {E m } m= o,..., 2n +i C M and assume the ordering (5.56). Then 
the isospectral manifold of smooth bounded real-valued solutions u G C 00 (R)nL 00 (R) 
of s-CH„(w) = is given by a real n-dimensional torus T n . 

Proof. The discussion in Remark 5.6 and Theorem 5.8, shows that the motion of 
each fi,j(x) on K, n proceeds topologically on a circle and is uniquely determined by 
the initial data fik(xo), k = 1, . . . , n. More precisely, the initial data 

p,j(x ) = {vj{x ),y{fj.j(x ))) = (fj, j (xo),-Hj(xo)G n (fj, j (xo),x )), ^ 6 ^ 

Hj{x ) e [E 2j -i,E 2 j], j = l,...,n 

are topologically equivalent to data of the type 

(/x j (xo),o- j (x )) G [E 2j -i,E 2j ] x {+,-}, j = l,...,n, (5.68) 

the sign of <Tj(xo) depending on /tj(xo) G II-|-. If some of the /itfc(xo) G {E 2 k-i, E 2 k}, 
then the determination of the sheet II± and hence the sign crfe(xo) in (5.68) becomes 
superfluous and is eliminated from (5.68). Indeed, since by (2.16), 

Aij(x ) 2 G„(Mj(xo),x ) = R 2n+2 (Hj( x o)), (5-69) 

G n (/j,j(xo), xo) is determined up to a sign unless Hj(xo) hits a spectral gap endpoint 
E 2 j_i, E 2 j in which case G„(^ :) (xo), xo) = i?2n+2(A*j(xo)) = and the sign ambigu- 
ity disappears. The n data in (5.68) (properly interpreted if Hj(xo) G {E 2 j_\, E 2 j}) 
can be identified with circles. Since the latter are independent of each other, the 
isospectral manifold of real- valued, smooth, and bounded solutions of s-CH„(u) = 
is given by a real n-dimensional torus T". □ 

Remark 5.10. 

(i) For simplicity we only focused on the case Au — u xx > 0. The opposite case 
Au — u xx < is completely analogous and results in a reflection of E m , to = 
0, . . . , 2n + 1, and Hj(x), Vj(x), j = 1, . . . ,n, about z = 0, etc. 

(ii) The time-dependent case also offers nothing new. Higher-order CH r flows drive 
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each fij(x,t r ) around the same circles as in the stationary case in complete analogy 
to the familiar KdV case. 

In summary, one observes that the reality problem for smooth bounded solutions 
of the CH hierarchy, assuming the ordering (5.56) (resp., the one obtained upon 
reflection with respect to z — 0), parallels that of the KdV hierarchy with the basic 
self-adjoint Lax operator (the one-dimensional Schrodinger operator) replaced by 
the self-adjoint Hamiltonian system (5.6). 

The following result was found in response to a query of Igor Krichever, who 
inquired about the significance of the eigenvalue ordering (5.56) (or the one obtained 
upon reflection at z = 0). As it turns out, such an ordering is crucial if one is 
interested in smooth algebro-geometric solutions on R. 

Theorem 5.11. Suppose 

E m < E m+1 , m = 0, . ..2n, 

and E 2 j -i = (resp., E 2 j = 0) for some j E {1, . . . , n}. (5.70) 
Fix Xq E K, and assume that the initial data 

{fij{x ) = {/J,j(x ),-iJ,j(x )G n (nj(x ),x ))}j= lt ... tn C K, n (5.71) 
for the Dubrovin equations (3.31) are constrained by 
Hj(x ) E [E 2 j-i,E 2 j], j e {1, . . . ,n}\{j }, 

and /J. jo {x ) G {E 2jo -!, E 2jo ] (resp., n jo (x ) E [E 2jo -u E 2jo )). (5.72) 
Then there exists a set fl^ C M of the type 

M = K\{£ fe } fceZ , £ fe < £ fe+1 , k G Z, lim £ fe = -oo, lim £ fe = oo, (5.73) 

k[ — oo k^oo 

such that the Dubrovin initial value problem (3.31), (5.57), (5.58) has a unique 
solution {Aj}j=i,...,n C K, n satisfying 

fter^.q, i = l,...,n, (5.74) 

and i/ie projections [ij remain trapped in the intervals [E 2 j-\, E 2 j], j = l,...,n, 
for all x E O^, 

/ij(ar) G [S 2 j-i,^2j], j = 1, • • • ,n, x G ft^. (5.75) 
Moreover, u defined by the trace formula (3.34), that is, 

^ n _^ 2n-t-l 
j — 1 m— 

satisfies 

u E C°°(Q.^), u is real-valued, (5-77) 

4u - u XI > 0, a; G ty,, (5.78) 

and </ie ni/i stationary CH equation 

s-CH„(it) = on n M , (5.79) 

wifft integration constants C£ in (5.65) given by a = ce(E), I = 1, . . . ,n, according 
to (2.26), (2.27). At each u x exhibits a singularity of the type 

u x (x) = C k (x - ik)- 1 '* + o((x-Ck)- 1/3 ), C fc ^0, fcGZ. (5.80) 
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In particular, u £ C 1 (R), u x ^ L°°(E). T/ie isospectral manifold corresponding to 
(5.70)-(5.72) is then given by T™ _1 x E. 

Proof. One can follow the proof of Theorem 5.8 and Corollary 5.9 with one im- 
portant twist, though, since the right-hand side of the Dubrovin equation of /j,j 
blows up as Hj — > E 2 j -i (rcsp., E 2 j ). For notational convenience and without 
loss of generality we may assume E\ — (and hence jo = 1) in the following. 
Recalling the Dubrovin equations (3.31), one verifies that its solutions are smooth 
with respect to x as long as \i\ stays away from E\ = 0. Varying x G E, the sign 
restrictions on fi\ tX in terms of the right-hand side of the corresponding equation 
in (3.31) eventually accelerate fi\ into E\ = as x tends to some G E, and we 
now analyze what happens to all fij for x in a neighborhood of Recalling the 
local coordinate crz 1 / 2 , a = ±1, near E\ = and hence introducing 

Ci(aO = cr(— /^(a;)) 1 / 2 for Mi(^) sufficiently close to E\ = as x — > (5.81) 

and the corresponding point /ti = — £ 2 G /C n , the Dubrovin equation for becomes 
for x near £ fe , 



£,*(*) = gyg^ ; II(-Ci( a; ) 2 - WW)" 1 , (5-82) 
CiCi(x)- 2 (l + o(l)), 



fe (x) = 2^n( ft (,)- w ( I ))- 1 , i = 2,...,n. (5.83) 

for some constant C\ ^ 0. (Here we implicitly assume that no other fij, j = 2, . . . , n 
simultaneously hits E 2 j-i or E 2 j as x — > Otherwise one simply resorts to the 
proper local coordinate for such a fij. We omit the details.) To treat the singularity 
°f Ci,x as (i(x) — > for a; — > wc now resort to a well-known trick described, 
for instance, in [15, Theorem 3.2.2] in the context of scalar first-order differential 
equations. Instead of looking for solutions Ci, fij as functions of x, we now look 
for x — x(d), jlj = /5j(Ci) as functions of Ci, where we denote Mj(Ci) = f-j{x), 
j = 2,...,n. Then (5.82), (5.83) turn into 

z'(Ci) = -4L- I](-Ci 2 - A/(Ci)), (5-84) 

y(-C{) 1=2 

a jtCl (x) = 2 y -^&Jf- n(-Ci 2 - a/(Ci)) n^(Ci) - MCi)r\ (5.85) 

eft 

j = 2, . . . ,n. 

Since the right-hand sides in (5.84), (5.85) are holomorphic with respect to the 
n variables CbM2, ■ ■ ■ ,Mn f° r Ci nea r zero and jlj near [E 2 j_i, E 2 j], j = 2, . . . ,n, 
equations (5.84) and (5.85) yield solutions x, ji 2 , . . . , jl n holomorphic with respect 
to £1 near £1 = 0. In particular, since 

*'(Ci) = Ci/^+OiCf), MKCi) = n ^Ci 2 +0(Ci 3 ), J = 2,...,n, (5.86) 
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for some constants Cj 7^ 0, j = 1, . . . , n, one obtains 

x(Ci) = Jk + Cf/i^+OiCt), (5.87) 

fc(Ci) = Ai(0) + QCi 3 /3 + 0(Ci 4 ) = + QCi 3 /3 + 0(Ci 4 )- (5-88) 
Ci— >o 

Thus, inverting #((1), one observes 

Ci(x) = [3d(x-6)] 1/3 + 0((x-6) 2/3 ), (5.89) 

a;— >5fc 

/^(z) = ^(a)+C 1 C ] (x-a) + 0((x-Ck) 4/3 ), j = 2,...,n, (5.90) 
and hence, 

= -(2/3)(3C 1 ) 2 / 3 (x-a)- 1/3 +o((x-a)- 1/3 ) (5.91) 

and 

n 

u x (x) = {1/2)^2 n jiX (x) 
i=i 

= -(l/3)(3C 1 ) 2 / 3 (x-a)- 1/3 + o(( 2 ;-a)- 1/3 ). (5.92) 

The singular behavior (5.91), (5.92) repeats itself after each revolution of [i\ around 
its circle and occurs whenever /zi passes again through E\ = 0, giving rise to 
the exceptional set {Cfcjfcez in (5.73). Hence, fij £ C 1 (M), j = 2, ...,n, while 
jii /x = — Ci(a^) 2 blows up whenever x approaches an element of {£fc}fcez- The rest of 
the discussion follows as in Theorem 5.8 and Corollary 5.9. Since fii(xo) = E\ = 
is not an admissible initial condition in (5.72), one point must be removed from the 
circle associated to /ii which topologically results in R instead of S 1 and hence in 
the noncompact isospectral manifold T"^ 1 xR. □ 

Thus, smooth algebro-geometric CH solutions require Eq = or i?2n+i = 0. 

Acknowledgments. We are indebted to Mark Alber, Roberto Camassa, Adrian 
Constantin, Yuri Fcdorov, and Igor Krichever for stimulating discussions. 
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